Grid simulation services for the medical community

The first part of this paper presents a selection of medical simulation applications, including image reconstruction, near real-time registration for neuro-surgery, enhanced dose distribution calculation for radio-therapy, inhaled drug delivery prediction, plastic surgery planning and cardio-vascular system simulation. The latter two topics are discussed in some detail. In the second part, we show how such services can be made available to the clinical practitioner using Grid technology. We discuss the developments and experience made during the EU project GEMSS, which provides reliable, efficient, secure and lawful medical Grid services

The Epitheliome: agent-based modelling of the social behaviour of cells

We have developed a new computational modelling paradigm for predicting the emergent behaviour resulting from the interaction of cells in epithelial tissue. As proof-of-concept, an agent-based model, in which there is a one-to-one correspondence between biological cells and software agents, has been coupled to a simple physical model. Behaviour of the computational model is compared with the growth characteristics of epithelial cells in monolayer culture, using growth media with low and physiological calcium concentrations. Results show a qualitative fit between the growth characteristics produced by the simulation and the in vitro cell models

Virtual Coronary Intervention: A Treatment Planning Tool Based Upon the Angiogram

Objectives: This study sought to assess the ability of a novel virtual coronary intervention (VCI) tool based on invasive angiography to predict the patient's physiological response to stenting. Background: Fractional flow reserve (FFR)-guided percutaneous coronary intervention (PCI) is associated with improved clinical and economic outcomes compared with angiographic guidance alone. Virtual (v)FFR can be calculated based upon a 3-dimensional (3D) reconstruction of the coronary anatomy from the angiogram, using computational fluid dynamics (CFD) modeling. This technology can be used to perform virtual stenting, with a predicted post-PCI FFR, and the prospect of optimized treatment planning. Methods: Patients undergoing elective PCI had pressure-wire-based FFR measurements pre- and post-PCI. A 3D reconstruction of the diseased artery was generated from the angiogram and imported into the VIRTUheart workflow, without the need for any invasive physiological measurements. VCI was performed using a radius correction tool replicating the dimensions of the stent deployed during PCI. Virtual FFR (vFFR) was calculated pre- and post-VCI, using CFD analysis. vFFR pre- and post-VCI were compared with measured (m)FFR pre- and post-PCI, respectively. Results: Fifty-four patients and 59 vessels underwent PCI. The mFFR and vFFR pre-PCI were 0.66 ± 0.14 and 0.68 ± 0.13, respectively. Pre-PCI vFFR deviated from mFFR by ±0.05 (mean Δ = -0.02; SD = 0.07). The mean mFFR and vFFR post-PCI/VCI were 0.90 ± 0.05 and 0.92 ± 0.05, respectively. Post-VCI vFFR deviated from post-PCI mFFR by ±0.02 (mean Δ = -0.01; SD = 0.03). Mean CFD processing time was 95 s per case. Conclusions: The authors have developed a novel VCI tool, based upon the angiogram, that predicts the physiological response to stenting with a high degree of accuracy

Effect of side branch flow upon physiological indices in coronary artery disease

Recent efforts have demonstrated the ability of computational models to predict fractional flow reserve from coronary artery imaging without the need for invasive instrumentation. However, these models include only larger coronary arteries as smaller side branches cannot be resolved and are therefore neglected. The goal of this study was to evaluate the impact of neglecting the flow to these side branches when computing angiography-derived fractional flow reserve (vFFR) and indices of volumetric coronary artery blood flow. To compensate for the flow to side branches, a leakage function based upon vessel taper (Murray’s Law) was added to a previously developed computational model of coronary blood flow. The augmented model with a leakage function (1Dleaky) and the original model (1D) were then applied to predict FFR as well as inlet and outlet flow in 146 arteries from 80 patients who underwent invasive coronary angiography and FFR measurement. The results show that the leakage function did not significantly change the vFFR but did significantly impact the estimated volumetric flow rate and predicted coronary flow reserve. As both procedures achieved similar predictive accuracy of vFFR despite large differences in coronary blood flow, these results suggest careful consideration of the application of this index for quantitatively assessing flow

Estimation of valvular resistance of segmented aortic valves using computational fluid dynamics

Aortic valve stenosis is associated with an elevated left ventricular pressure and transaortic pressure drop. Clinicians routinely use Doppler ultrasound to quantify aortic valve stenosis severity by estimating this pressure drop from blood velocity. However, this method approximates the peak pressure drop, and is unable to quantify the partial pressure recovery distal to the valve. As pressure drops are flow dependent, it remains difficult to assess the true significance of a stenosis for low-flow low-gradient patients. Recent advances in segmentation techniques enable patient-specific Computational Fluid Dynamics (CFD) simulations of flow through the aortic valve. In this work a simulation framework is presented and used to analyze data of 18 patients. The ventricle and valve are reconstructed from 4D Computed Tomography imaging data. Ventricular motion is extracted from the medical images and used to model ventricular contraction and corresponding blood flow through the valve. Simplifications of the framework are assessed by introducing two simplified CFD models: a truncated time-dependent and a steady-state model. Model simplifications are justified for cases where the simulated pressure drop is above 10 mmHg. Furthermore, we propose a valve resistance index to quantify stenosis severity from simulation results. This index is compared to established metrics for clinical decision making, i.e. blood velocity and valve area. It is found that velocity measurements alone do not adequately reflect stenosis severity. This work demonstrates that combining 4D imaging data and CFD has the potential to provide a physiologically relevant diagnostic metric to quantify aortic valve stenosis severity

Adaptation and development of software simulation methodologies for cardiovascular engineering: present and future challenges from an end-user perspective

This paper describes the use of diverse software tools in cardiovascular applications. These tools were primarily developed in the field of engineering and the applications presented push the boundaries of the software to address events related to venous and arterial valve closure, exploration of dynamic boundary conditions or the inclusion of multi-scale boundary conditions from protein to organ levels. The future of cardiovascular research and the challenges that modellers and clinicians face from validation to clinical uptake are discussed from an end-user perspective

MRI model-based non-invasive differential diagnosis in pulmonary hypertension

Pulmonary hypertension(PH) is a disorder characterised by increased mean pulmonary arterial pressure. Currently, the diagnosis of PH relies upon measurements taken during invasive right heart catheterisation (RHC). This paper describes a process to derive diagnostic parameters using only non-invasive methods based upon MRI imaging alone.Simultaneous measurements of main pulmonary artery (MPA) anatomy and flow are interpreted by 0D and 1D mathematical models, in order to infer the physiological status of the pulmonary circulation. Results are reported for 35 subjects, 27 of whom were patients clinically investigated for PH and eight of whom were healthy volunteers. The patients were divided into 3 sub-groups according to the severity of the disease state, one of which represented a negative diagnosis (NoPH), depending on the results of the clinical investigation, which included RHC and complementary MR imaging.Diagnostic indices are derived from two independent mathematical models, one based on the 1D wave equation and one based on an RCR Windkessel model. Using the first model it is shown that there is an increase in the ratio of the power in the reflected wave to that in the incident wave (Wpb/Wptotal) according to the classification of the disease state. Similarly, the second model shows an increase in the distal resistance with the disease status. The results of this pilot study demonstrate that there are statistically significant differences in the parameters derived from the proposed models depending on disease status, and thus suggest the potential for development of a non-invasive, image-based diagnostic test for pulmonary hypertension

CDF- and Bernoulli-based pressure drop estimates: A comparison using patient anatomies from heart and aortic valve segmentation of CT images

Purpose: An aortic valve stenosis is an abnormal narrowing of the aortic valve (AV). It impedes blood flow and is often quantified by the geometric orifice area of the AV (AVA) and the pressure drop (PD). Using the Bernoulli equation, a relation between the PD and the effective orifice area (EOA) represented by the area of the vena contracta (VC) downstream of the AV can be derived. We investigate the relation between the AVA and the EOA using patient anatomies derived from cardiac computed tomography (CT) angiography images and computational fluid dynamic (CFD) simulations. Methods: We developed a shape-constrained deformable model for segmenting the AV, the ascending aorta (AA) and the left ventricle (LV) in cardiac CT images. In particular, we designed a structured AV mesh model, trained the model on CT scans, and integrated it with an available model for heart segmentation. The planimetric AVA was determined from the cross-sectional slice with minimum AV opening area. In addition, the AVA was determined as the non-obstructed area along the AV axis by projecting the AV leaflet rims on a plane perpendicular to the AV axis. The flow rate was derived from the LV volume change. Steady-state CFD simulations were performed on the patient anatomies resulting from segmentation. Results: Heart and valve segmentation was used to retrospectively analyze 22 cardiac CT angiography image sequences of patients with non-calcified and (partially) severely calcified tricuspid AVs. Resulting AVAs were in the range of 1-4.5cm2 and ejection fractions (EFs) between 20-75%. AVA values computed by projection were smaller than those computed by planimetry, and both were strongly correlated (R2 = 0.995). EOA values computed via the Bernoulli equation from CFD-based PD results were strongly correlated to both AVA values (R2 = 0.97). EOA values were ˜ 10% smaller than planimetric AVA values. For EOA values < 2.0cm2, the EOA was up to ˜ 15% larger than the projected AVA. Conclusions: The presented segmentation algorithm allowed to construct detailed AV models for 22 patient cases. Because of the crown-like 3D structure of the AV, the planimetric AVA is larger than the projected AVA formed by the free edges of the AV leaflets. The AVA formed by the free edges of the AV leaflets was smaller than the EOA for EOA values < 2.0cm2. This contradiction with respect to previous studies that reported the EOA to be always smaller or equal to the geometric AVA is explained by the more detailed AV models used within this study